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ABSTRACT 

A mean-field model is proposed as a test case for tricritical series analyses meth- 
ods. Derivation of the 50**^ order series for the magnetization is reported. As the first 
application this series is analyzed by the traditional slicewise Pade approximant method 
popular in earlier studies of tricriticality. 
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1. INTRODUCTION 



Development of algorithms for the numerical investigation of tricritical behavior in 
a two-variable phase diagram has been an elusive goal for many years in the context 
of both magnetism and polymer studies. Two such systems that have been the focus 
of special interest are random field models for which the tricritical point is expected 
in the two-parameter space of temperature and randomness strength [1], and the 9- 
point transition of linear polymers [2] for which the second variable is the "stickiness" 
fugacity leading to collapse. For these and other systems with tricritical points, all the 
standard large-scale numerical methods have been utilized: Monte Carlo [3,4,5] and 
transfer matrix techniques [5,6] (the latter for 2D models), and series analyses [7]. 

There are many complicating factors for such studies. The models involved are 
rarely simple and must be analyzed at multiple points in the two-parameter space. A 
"test problem" with an exact solution and a series expansion (such as the 2D Ising 
model which is widely used in calibrating new techniques for second order transitions) 
has been missing until now for two-variable problems with tricritical points. 

Some of the technical problems are specific to the numerical technique. In the 
random field and ^-point transition models reaching true equilibrium in Monte Carlo 
simulations, for instance, is complicated by slowdown effects in the low-temperature or 
dense (collapsed) regions of the phase diagram. We shall not address simulations further 
in this paper but note that there is considerable controversy over tricritical behavior 
between different recent studies of the 2D collapse transition [4]: 9 vs. 9' points, etc. 

Series expansions do not suffer from equilibration problems, and in many cases 
lattice models are amenable to the generation of series for all points in the two- variable 
phase diagram. In particular, two- variable series have been developed for both Ising 
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random field (IS*'^ order in general dimension) [1] and several polymer problems [2]. 
However for series the main complication is that the very nature of the tricritical be- 
havior means that techniques for studying both first and second order transitions must 
be applied. While excellent techniques for identifying first order transitions via simula- 
tion have recently been developed [8], methods to identify first order transitions from 
series expansions are unreliable [9] unless both low and high temperature series exist. In 
one notable case (the FCC-lattice Blume-Capel model) where such expansions on both 
sides of the transition were developed [10] a satisfactory characterization of tricriticality 
was made from series. 

Despite the absence of a systematic, well tested, approach to studying tricriticality 
when series from only one temperature direction are available, some attempts to do 
so have been made. The variety of the makeup methods used, frequently resulted in 
differences in answers to the basic question concerning the existence of a tricritical point, 
not to speak of exponent and other parameter estimates which could be attributed to 
the diversity of methods rather then to the quality of the series expansions available. 

This work reports two developments towards systematizing tricritical series analy- 
ses. Firstly, we derive a test series based on a mean-field model with a tricritical behavior 
which is well understood and has most of the features of "real" tricritical points of 2D 
and 3D systems. Secondly, we apply the standard "slicewise" Pade method to this new 
series. We identify those features of the Fade approximant approach which can be re- 
garded as signatures of a tricritical point in the phase diagram and which were noted 
in some early studies of tricriticality by series [7] . 

We conclude, however, that this most straightforward Pade method is not suit- 
able as an accurate and systematic general analyses technique, and it can be only used 
for exploratory studies or in conjunction with other information. Application of the 
slicewise Pade method to certain random-field model series will be reported in a forth- 
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coming publication. However, the door is still wide open for developing a systematic 
series-analysis method, possibly based on elaborations of the two-variable differential 
approximant techniques used successfully for bicritical points [11]. 

One interesting aspect of our test series derivation, reported in Section 2, and its 
Fade analysis in Section 4, is that both rely heavily on novel large-scale computational 
abilities. Series derivation required a large Mathematica run, whereas Fade analysis 
employing simultaneously many Fade-pole calculations and extensive graphics repre- 
sentations of the data, revealed new features not accessible to earlier studies from the 
seventies [7]. Thus, future tricritical series analyses are likely to be large-scale compu- 
tational projects. Section 3 summarizes the tricritical phase diagram of the mean-field 
model used in test-series studies. Finally, Section 5 is devoted to some concluding 
remarks and to acknowledgments. 
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2. THE MEAN-FIELD MODEL 



In this section we report derivation of a low-temperature two-variable series for a 
tricritical point in an infinite range model with mean-field critical behavior. There are 
two kinds of solvable Ising-type infinite-range models. In the first and more familiar 
type the spins are ±1, but their interaction energy which is the function of the total 
magnetization, is essentially arbitrary. Usually it is selected (or Taylor-expanded) as a 
polynomial in the total magnetization so that the resulting constrained free energy, as 
function of the magnetization, m, resembles the Landau expansion. 

The second type of a mean-field model [12] is defined by having a simple quadratic 
energy but a complicated entropy-like contribution due to assigning essentially arbitrary 
measure in the evaluation of the partition function of a system of scalar spins which 
vary in (— oo, -|-oo). 

Both types of infinite-range models suffer from the difficulty that the low-temperature, 
T, behavior is different from the short-range lattice models. Indeed, for short-range 
models the low-T series expansions are in terms of the Boltzmann factors of excitations 
above a reference ground state, of the form exp(— A£'/A;T), where = 0(1) is the 
energy cost due to a local structure (overturned spins, broken bonds, etc.). However, 
this "locality" of the excitation structure is lost for infinite-range models. The ±l-spin 
models have entropic contribution to the constrained (fixed-m) free energy with singu- 
larities of the type ~ (1 — m) ln(l — m) near magnetization m = 1 (and similarly near 
m = —1). Thus setting up a low-T expansion presents a mathematical challenge. 

The models of [12] are less troublesome in this respect. One can get a well-controlled 
series in powers of T itself. This "soft" T-dependence is an artifact of the infinite-range 
model, and there are some other artificial features near T = 0, but the series is well 
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defined and can be derived in closed form to any fixed order given sufficiently powerful 
computational facilities. Thus, we choose to work with the model of [12] here. 

The energy of the interacting scalar spins, CTj, is taken as 



E = 



J 

'2N 



N 



(2.1) 



where N is the number of spins, and J > 0. The partition function is defined as 



N 

, / exp{-E/kT)l[di^{ai) , 
i=i 



(2.2) 



where the spins are weighed with measure dfj,{a). The order parameter is obtained from 



N 



aieM-E/kT)Y[dii{ai) 



(2.3) 



i=l 



The Gaussian-integral method [12] can be used to show that in the limit A?" — > oo 
the free energy, f, in Z — exp{—Nf), can be obtained as 



/ 



mm 

X 



kTx^ 
2J 



-Q{x) 



(2.4) 



where 



Q{x) = lnj e'^diiia) . 



(2.5) 



If the minimum in (2.4) is at some x = Xm, then one can further show that 



m = 



dQ 
dx 



J 



(2.6) 



where the last equahty follows from the fact that the global minimum is obtained at 
one (or more) roots of 



dQ _ kTx 
dx~ J ■ 

Thus, we note that m = kTxm/J, i-e., m oc T for low temperatures. This is one of 
those artificial infinite-range model features. It turns out convenient to work with Xm 
directly rather than with to, as the order-parameter-like quantity for series analysis. Of 
course, the actual critical-tricritical-first-order behavior is at T > so the difference 
only affects the form of analytic corrections to scaling. 

In order to have a solvable model with tricritical behavior, we take Q{x) as an even, 
six-degree polynomial in x. There is still freedom in selecting the coefficients, etc. We 
choose to work with dimensionless parameters which, disregarding various dimensional 
factors, amounts to effectively putting 

J = k/2 , (2.8) 



Q{x) =x^ + {U -l)x^ ~x^' , (2.9) 
so that our choice corresponds to 

f = mm\(T -l)x'^ -(U -l)x'^ + x^] . (2.10) 

X 

This is the simplest Landau-expanded form to yield the tricritical point, at {T,U) = 
(1,1), in the two-parameter space of the (dimensionless) temperature, T, and another 
(dimensionless) "coupling constant," U. 
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On the low-T side there are two symmetric roots of (2.7), and there is always one 
root at a; = 0. We consider the root 

^ 0) th.e actual series is conveniently generated 

for 



iVs=\/ y/{ir- 1)2 - 3(T - 1) + C/ - 1 . (2.11) 



By utilizing Mathematica, we derived the order 50 double series in T and U for this 
order-parameter quantity. This series, i.e., the first 2601 coefficients c^, for i,j = 
0, . . . , 50, in 

V3x^ = 5^^CyrC/^' , (2.12) 

i=0 j=0 

can be obtained via electronic mail, from the authors, on request. 
We also derived the functions ai(U) in 

^/Zx^ = Y,(^^{U)T' , (2.13) 

i=0 

for z = 0, 1, . . . , 50. These functions are available in the FORTRAN form, via electronic 
mail. 
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3. THE TRICRITICAL PHASE DIAGRAM 



For U < 1, there is a second-order transition line at T = 1, at which the order 



portionahty constant diverges as — U for t/ — > 1 . At the tricritical point, the 

order-parameter vanishes according to ~ (1 — T)^/^. 

For U > 1 there is a first-order transition. The fine T = 1 stiU have special signifi- 
cance as the mean-field spinodal (see further below) above which the high-temperature, 
zero-order-parameter phase exists. However, it is not seen in the low-T expansion. The 
actual first-order transition line is determined by the condition that the minima at 
Xm > and at x = are equal; cf. (2.10). A somewhat lengthy calculation yields 



for the first-order transition line at U > 1. 

Along this line the low-T-side order parameter vanishes according to oc — 1 as 
t/ ^ 1+, i.e., on approach to the tricritical point. However, for fixed U > 1, the order 
parameter is finite at the first-order transition. 



from the low-T side, and it vanishes from the high-T side. 

In short-range Ising-type lattice models, if one attempts to analytically continue 
the thermodynamic functions "through" the first-order line, one encounters an essential 
singularity at the first-order transition. This singularity is due to droplet excitations; it 



parameter approaches zero according to the mean-field law Xm oc Vl — T. The pro- 




(3.1) 




(3.2) 
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is weak and its detection in series analysis has rarely been accomplished unambiguously 
[9] . Specifically, its manifestation within the traditional Pade method aimed at detecting 
power-law divergences, is at best indirect via sequences of weak, alternating poles and 
zeroes of the approximants [9]; see the following sections for further discussion. While 
the incorporation of this essential singularity must be ultimately a goal for a fully 
systematic series-analysis method of tricriticality, at the present state of the art and 
available series lengths, its presence will have little efi^ect in any series study. 

There is no essential singularity for infinite-range models as there are no droplet 
excitations, only uniform ones. When the thermodynamic quantities are continued 
past the first-order transition, one encounters a spinodal line at which the low-T order 
parameter Xm > ceases to be a local minimum of the constrained free energy. (While 
at the first-order line it ceases to be the global minimum.) This mean-field spinodal line 
is at 



The existence of a sharp spinodal-type singularity is an artifact of the infinite-range 
model. However, for short-range models traces of spinodal-type behavior have been 
noted in available-length series analyses [9]. These are artifacts of employing approxi- 
mants (within Pade or other analysis methods) which fit the data to a form suggestive 
of a sharp singularity; see [9] for further discussion. 

Near the tricritical point, one can write the low-T-side scaling form in terms of the 
scaling variables 



T=l + l(t/-l)2 . 



(3.3) 



t = T- 1< 



and 



u = U-l 



(3.4) 
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V3xmC^{-ty/'F_(u{-t)-'/^) , (3.5) 
where the scahng form (3.5) apphes for t, u — > and the scahng function is 

F_(C) = i/c + • (3.6) 
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4. SLICEWISE FADE ANALYSIS 



The slicewise Pade analysis is perhaps the simplest, single- variable approach to 
analyzing double-series expansions [7]. Thus, we calculate, for fixed U, approximants 
to the series coefficients in 

v^x^(T,C/-fixed) = ^a,(t/)TV (4.1) 
The coefficients Oj are approximated by 



aiC^Yl ^^i^^' ' (4.2) 



where in our case jmax = 50. 



The order [M/ L] dlog-Pade approximant to the derivative = dxm/ dT is defined 
as the rational approximant of the form 



^ Po+PiT + P2T^ + ...+PmT'' 3. 
x'^ l + qiT + 92^2 + . . . + qlT^ ■ ^ ■ ^ 

where the derivative rather than used in order to have both the numerator 

and denominator of the left-hand side diverge at their first singularity as T is increased 
from zero. Specifically, for fixed U, we have 



x'^^[TiU)-T]-^ , (4.4) 

with T(U) = Tc = 1, S = 1/2 at the critical line for U < 1, and T(l) = 1,5 = 3/4 at 
the tricritical value U — 1. For the infinite-range model T{U) equals the spinodal value 



(3.3), while B = 1/2, for fixed U > 1. Note that ioi U < 1 the exponent B is related to 
the order-parameter exponent usually denoted by /3, via B = 1 — j3. 

The coeflBcients po,...,M and qi,...,M are calculated in a standard fashion [13] to have 
the power series of the right-hand side of (4.3) reproduce the first M + L + 1 power-series 
coefficients of the left-hand side. Of special interest are the poles of the approximant 
which here depend parametrically on U and will be loosely denoted simply by T{U). 
It is anticipated that for power-law singularities (4.4) a "stable" pole location will be 
found in the highest-order, near diagonal approximants (i.e., M ~ L and M + L -\- 1 
close or equal to the order of the available series for x'^^/x'^) such that for T near T{U) 
the right-hand side of (4.3) approximates the behavior suggested by (4.4), 



x'^^T{U)-T ■ ^ • ^ 

Thus, the residue at the stable pole approximates the exponent in (4.4), —B. 

For essential singularities and other singularities associated with branch cuts, it has 
been noted [9] that Fade approximants sometimes yield a sequence of alternating weak 
poles and zeros (zeros of the denominator and numerator, respectively) which mimic the 
branch cut. However, the Fade method is well suited only for single- variable expansions 
of functions with power-law singularities. It is worth pointing out that recently exact 
results were derived for certain models of partially convex lattice vesicles [14] which 
show tricritical behavior with essential singularities at the "first-order line" of their 
phase diagram. Double- variable expansions can be derived for these models [15] and 
possibly used as test series for methods to detect essential singularities, etc. However, 
we note that the singularities of the models of [14,15] seem to be natural-boundary- 
type and differ from droplet-type essential singularities anticipated at Ising first-order 
transitions. 



It is also important to point out that in the shcewise Pade method used here 
the approximation if two-fold: the coefficients Oj are calculated approximately via the 
truncated series (4.2); the Pade method is applied to the truncated series (4.1). We 
kept the order of the Pade approximation, M + L + 1, at about half the order of the 
truncation in (4.2), which is jmax = 50. Still, as examples below illustrate, the method 
fails for U^l which is presumably due to the truncation (4.2). While the truncation 
(4.2) is the approach used in the early literature, it is natural to consider improvement 
of this approximation: we address this issue later in this section. 

As our first example, it is useful to consider approximants with only a single pole. 
For instance. Figure 1 shows the pole location for the approximant [19/1]. We note 
that the approximant becomes "defective" for U^l in that the location of the pole 
has nothing to do with the actual spinodal line. Figure 2 also illustrates that the 
residue provides a poor approximation to the exponent B. The latter, however, can 
be blamed on our use of the extremely off-diagonal approximant. This "hooking" of 
the approximant away from the actual phase-transition line has been noted in earlier 
studies [7]. The hooking can also be in the direction opposite to that of Figure 1, as 
occurs, for instance, in the [30/1] approximant not shown here. 

Study of diagonal and also numerous near-diagonal approximants (only two are ac- 
tually illustrated here; see Figures 3, 4, 5) reveals that they indeed significantly improve 
the exponent B estimate for small U < 1. Most approximants also yield B values quite 
close to 3/4 at U = 1. However, the phase diagram is not well-represented near and 
above the tricritical point. In the wide crossover regime near U = 1, the approximants 
become defective: the "physical" pole alternates among several branches of the roots 
of the denominator; see Figure 3 for the [10/10] approximant. While these branches 
first roughly follow the exact spinodal location above U = 1, they soon "hook away." 
The exponent estimates also show a rather irregular crossover pattern near U = 1; 
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see for instance Figure 4 for [10/10]. The approximation fully deteriorates soon above 
U = 1. Furthermore, it seems that the quality of the approximation is affected little by 
increasing the approximant order as illustrated by the case [17/17] in Figure 5. 

Our conclusion is that the slicewise Pade method as implemented [7], can at best 
provide a qualitative indication of the presence of a tricritical point. A wide region of 
irregular, "defective" approximant behavior develops on approach to tricriticality. The 
hooking noted in earlier studies [7] can be attributed to changeovers among the various 
pole branches. Thus plotting all the poles and all the residues allows a rough location 
of the tricritical point and estimation of the exponent. But otherwise the slicewise 
Pade method should not be regarded as a systematic technique for analyzing tricritical 
behavior. 

We now turn to the approximation involved in the truncation (4.2). Some series 
actually have the functions ai(U) as polynomials in U (or other appropriate expansion 
parameter). However, generally the double-series are available in the form truncated in 
both variables. It is natural to assume that the quality of the overall approximation 
can be improved by using resummation methods for the single- variable series (4.2). To 
our knowledge, no systematic procedures were developed in the literature; one could 
contemplate Pade or other resummation methods. However, in this work we limited 
ourselves to the following: we repeated the preceding analysis with the exact functions 
ai{U), derivation of which was described in Section 2. 

Figure 6 illustrates the behavior of the poles for the case of the [10/10] Pade ap- 
proximant. It should be compared with Figure 3. We note that without the truncation 
(4.2), the approximation for U^l improves. Specifically, the leading pole now clearly 
follows the spinodal line for U > 1. However, the main problems remain: the crossover 
region near U = 1 is marked by the "defective" behavior where poles approach each 
other and switch role. Furthermore, the leading pole for C/ > 1 is accompanied by a 
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sequence of weak poles which mimic branch-cut effect. Figure 7 shows the residues at 
the poles of the [10/10] Pade approximant. It should be compared with Figure 4. The 
exponent estimation is only accurate for C/ < 1. In the wide crossover region near U = 1 
the defective behavior spoils the accuracy of the approximation, while for t/ > 1, the 
quality of the approximation is low presumably due to the accompanying weak poles. 
A similar behavior was observed for other near-diagonal Pade approximants calculated 
with the exact ai{U) values, and as before there was no visible improvement when the 
order, [M/ L] , of the approximant was increased. 
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5. CONCLUDING REMARKS 



As discussed in the previous section, the simplest Pade approach fails in several 
aspects near tricritical points. Let us now consider a "wish list" for a more systematic 
series analysis method. Firstly, we would expect to produce a smoothed out but regular 
approximation to the exponent, i.e., the spiked line in Figure 2, and to the phase 
diagram, i.e., to the phase-transition lines shown in Figure 1. The approximants should 
be sharpening up with the increased order of approximation. Secondly, we would also 
like to estimate the scaling form at the tricritical point, cf. (3.5), as in the differential- 
approximant analyses of bicritical points [11]. 

An added complication at tricritical points is the presence of singularities at the 
first-order transition line, as well as possible pseudo-singular spinodal behavior. Specif- 
ically, for Ising-type models there are weak essential singularities. However, for other 
models with soft-mode excitations, power-law spin-wave-type singularities are present 
at the first-order transition. For yet another class of models, the droplet picture may 
not be fully understood, such as for certain systems with randomness, or the first-order- 
regime singularities may not have been carefully discussed in the literature, such as for 
the polymer collapse. 

We note that the deterioration of the approximant quality in the slicewise Fade 
method has always occurred at larger U values. Improvement accomplished by avoiding 
the approximate truncation (4.2) was not sufficient for really accurate results. Other 
approaches involve, for instance, "slicing" along curvilinear paths that originate at the 
origin of the (T, U) plane, etc. These possibilities will be explored in future publications. 

All the above remarks indicate that series analysis of tricritical behavior promises 
to become an interesting and active field with the availability of new, long series and 
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modern computational facilities. It is hoped that the approach presented here will yield 
useful test series for these studies. 
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FIGURE CAPTIONS 



Figure 1: The symbols o, merging for most part as heavy sohd lines, show the pole 
T{U) of the [19/1] dlog-Pade approximant to the T-derivative of the magnetization 
series. The thin solid lines correspond to the first-order transition (lower curve, 
equation (3.1)), and to the mean-field spinodal line (upper curve, equation (3.3)), 
ioi U > 1. The tricritical point is at T = 1, C/ = 1, while for t/ < 1 there is the 
second-order transition at T = 1. 

Figure 2: The symbols o show the exponent B{U) as calculated from the residue at 
the pole of the [19/1] Pade approximant. The horizontal thin solid line corresponds 
to the exact value B — 1/2 ior U ^ 1. The vertical spike at U — 1 goes up to the 
exact value S(l) = 3/4. 

Figure 3: Shown are all the poles T(U) of the [10/10] Pade approximant which fit 
within the figure range. 

Figure 4: Exponent estimates B{U) from the residues of all the poles of the [10/10] 
Pade approximant. 

Figure 5: Shown are all the poles T{U) of the [17/17] Pade approximant. 

Figure 6: Shown are the poles T{U) of the [10/10] dlog-Pade approximant calculated 
with the exact coefficient values ai{U); cf. Figure 3. For U > 1, the thin solid 
line corresponds to the first-order transition (equation (3.1)), while the leading 
approximant (the smallest T{U) values) follows the mean-field spinodal (equation 
(3.3)), with the difference smaller than the size of the symbols. 

Figure 7: The exponent B{U) as calculated from the residue at the pole of the [10/10] 
Pade approximant with the exact ai{U) values; cf. Figure 4. The horizontal thin 
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solid line corresponds to the exact value B = 1/2 ior U ^1. The vertical spike at 
U =1 goes up to the exact value B{1) — 2,/A. 
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